The jamming transition and new percolation universality classes in particulate 

systems with attraction 
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We numerically study the jamming transition in particulate systems with attraction by investi- 
gating their mechanical response at zero temperature. We find three regimes of mechanical behavior 
separated by two critical transitions — connectivity and rigidity percolation. The transitions belong 
to different universality classes than their lattice counterparts, due to force balance constraints. We 
also find that these transitions are unchanged at low temperatures and resemble gelation transitions 
in experiments on colloidal and silica gels. 



The provocative conjecture that mechanical response 
at zero temperature is linked to slow dynamics at non- 
zero temperature in repulsive glassy systems has sparked 
tremendous interest in the jamming transition at "point 
J" [H. At point J athermal, frictionless systems with 
finite-range repulsive interactions [2[ undergo a single 
critical transition from an unjammed state with a van- 
ishing static shear modulus to a jammed state with a 
non-zero static shear modulus as density is increased. At 
the transition, it is observed that the average number of 
particles participating in the connected network jumps 
discontinuously from to L d (where d is the spatial di- 
mension and L is the system size), which makes the tran- 
sition first-order in the network order parameter. Recent 
work has indeed shown that properties of the jammed 
state are closely related to the slow dynamics in highly 
compressed repulsive glasses [3|- 

How does the jamming transition change for partic- 
ulate systems with attraction? It is likely that attrac- 
tive interactions will qualitatively change the nature of 
the jamming transition. For example, thermal systems 
with attraction can form repulsive glasses, but they can 
also form gels and attractive glasses [4]. Moreover, at 
zero temperature the mechanical properties of attractive 
granular materials and powders are quite different than 
those for dry granular media with purely repulsive inter- 
actions While the jamming phase diagram of Ref. [l[ 
needs revision to include attraction, the close correspon- 
dence between mechanical properties at zero tempera- 
ture and dynamics at non-zero temperature is robust. 
This has been demonstrated in experiments on attrac- 
tive colloidal suspensions |6[ where, as expected from the 
jamming phase diagram, gelation occurs upon increasing 
density, decreasing thermalization, and decreasing stress. 

In this Letter, we explore the jamming transition in 
attractive particulate systems by studying their mechan- 
ical response at zero temperature. A central conclusion 
from our work is that repulsive jamming is fundamentally 
different than attractive jamming, even for an infinitesi- 
mal amount of attraction. Instead of a single first-order 
transition in the purely repulsive systems, we observe 
two second-order transitions in the attractive systems — 



connectivity and rigidity percolation. These two transi- 
tions separate three distinct types of mechanical response 
and exhibit critical exponents that differ from those mea- 
sured in the connectivity and rigidity percolation transi- 
tions without force balance constraints. The transitions 
we observe at zero temperature are also present at small 
but finite temperatures and resemble gelation. 

Simulation procedure We investigate the quasistatic 
compression of an attractive particulate system in two 
and three dimensions at zero temperature. The system 
consists of N particles interacting via a pairwise, spher- 
ically symmetric potential, with a finite repulsive core 
and finite-range attraction. Simulations begin with a 
dilute collection of N spherical particles of diameter G{ 
randomly placed in a cubic cell with periodic boundary 
conditions. In each simulation step the diameters of all 
particles are increased by a small factor and then the 
potential energy is minimized using a conjugate gradient 
method. The force F{rij) between a pair of particles i 
and j depends on their separation relative to the sum 
of their radii = (pi +Oj)/2, and is plotted in the inset 
of Fig.pja). Since the overall scale Y of the force is irrel- 
evant at zero temperature, there are three independent 
parameters: the packing fraction </> (calculated using the 
location of the minimum in the force law), the range of 
the attraction and the minimum scaled force —C. 

Mechanical Response At zero temperature, whether 
a particulate system with attractions is jammed depends 
on two important factors: a) does the system contain 
a percolating cluster and b) if so, does the percolating 
cluster possess any floppy modes? A floppy mode Q is 
an infinitesimal deformation of the material that does 
not increase its energy. The number of floppy modes is 
equal to the number of non-trivial zero eigenvalues of the 
dynamical matrix for a given set of particle positions and 
radii [8]. If there are no non- trivial [9] zero eigenvalues, 
the potential energy increases with any deformation and 
all elastic constants are non-zero. We refer to this state 
as "jammed". If there is at least one non-trivial zero 
eigenvalue, the system is termed "unjammed". Note that 
this is a very strict definition of jamming. 

It is impossible to have a jammed system that does not 
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FIG. 1: (Color online), (a) The pressure P vs. (f) for two 
different compression procedures, which both begin at P — 
and cj) = 0.2. Red crosses correspond to the jamming path, 
where the system is compressed until it becomes jammed at 
the onset of rigidity percolation 0r. Black circles correspond 
to the connectivity percolation path, where the system is com- 
pressed to the onset of connectivity percolation 0p and then 
decompressed. The inset shows the central force law used in 
our simulations, (b) and (c) Screenshots at 0r for two dif- 
ferent systems. The blue-colored particles comprise the per- 
colating backbone, (d) Screenshot at 0p, with particles in 
the percolating cluster colored green, (e) Screenshot of the 
system pictured in (d) after decompression, where the perco- 
lating structure is stretched to its fullest extent. 

percolate, which leaves three distinct possibilities for its 
mechanical state: (1) nonpercolated, (2) percolated but 
unjammed, and (3) jammed. In repulsive systems only 
(1) and (3) occur and they are separated by the first- 
order jamming transition at point J. In the attractive 
particulate systems investigated here, all three mechani- 
cal states occur. A second-order connectivity percolation 
transition separates (1) and (2), and a second-order rigid- 
ity percolation transition separates (2) and (3). 

Properties of the three mechanical states are illustrated 



by measuring the pressure P. In Fig. [TJa) we plot P vs. 
(j) for two compression procedures in a d = 2 system with 
TV = 512, C = 10" 2 , and i = 0. In the first procedure, 
marked by red crosses, the system is compressed until 
it jams at 0r, after which the pressure becomes non- 
zero. At (/>r a percolating backbone forms with no floppy 
modes, also known as rigidity percolation [16]. Screen- 
shots of the system at </>r are shown in Figs. HJb) and 
(c). Note that while there are no floppy modes in the 
percolating backbone, there are trivial floppy modes in 
the system involving particles not in the backbone [l(j- 

In the second procedure, marked by black circles, the 
same initial system is compressed until a percolating net- 
work forms at (/>p, as illustrated in Fig. [lfd). Then the 
system is decompressed. We find that upon decompres- 
sion the pressure remains zero until a very low packing 
fraction. At this point the network is stretched to its 
fullest extent with no tension, as illustrated in Fig. HJe). 
Further decompression causes tension (P < 0) until the 
cluster becomes unstable, breaks, and P = 0. The forma- 
tion of a percolating cluster gives rise to novel mechanical 
behavior — infinitesimal deformations do not alter the po- 
tential energy but large deformations do. Thus, the me- 
chanical response of the percolated but unjammed state 
cannot be described by linear elasticity. 

The mechanical behavior described above is generic. 
For all compaction procedures we investigated, connec- 
tivity percolation and rigidity percolation are distinct 
transitions, with (j)p < (j)R. This occurs despite the fact 
that the mechanical response of the system is highly his- 
tory dependent. Since breaking a contact between parti- 
cles has an energy cost, contacts gained upon compres- 
sion are not necessarily lost upon decompression and we 
observe that the number of contacts per particle z never 
decreases when P = 0. Because contacts are not bro- 
ken, clusters of particles grow, or coarsen. This leads to 
work-hardening and lowers </>p and 0r. 

Connectivity Percolation Connectivity percolation 
occurs when a cluster of interacting particles first spans 
the system. The transition is characterized by three 
independent critical exponents, unless the hyperscaling 
assumption holds, in which case the number is two 
[ill ]. The choice of three independent exponents from a 
larger number of measurable exponents is arbitrary. We 
choose to focus on (1) the correlation length exponent v 
which quantifies the divergence of the linear extent £ of 
the largest non-percolating cluster £ oc \<\> — </>p | _zy , (2) 
the cluster-statistics exponent r which characterizes the 
number n s of clusters with s particles at the percolation 
transition n s ((j)p) oc s _r , and (3) the fractal dimension 
D at (j)p. The hyperscaling relation D(r — 1) = d holds 
if £ is the only relevant length scale in the system. 

The critical exponents of the connectivity percolation 
transition have been determined for lattices 11] and 
for continuum particulate systems without interparticle 
forces 12, [l3j]. The critical exponents, which are listed 
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FIG. 2: Measurement of the connectivity percolation expo- 
nents in d = 2. Circles (stars) correspond to monodisperse 
(bidisperse) distributions, (a) The square width of connec- 
tivity percolation thresholds A 2 as a function of system size 
N. The slope gives —2/dv. (b) The number of clusters of 
size s at <pp. The slope gives — r. (c) The average number of 
particles in the percolation cluster N per at 0p, as a function 
of system size N. The slope gives D/d. The best fit slopes 
determined from our analysis are included in each plot. 

in Table H are identical for lattice and continuum perco- 
lation and do not depend on the lattice type. In static 
particulate systems, force balance is crucial and has not 
been included in previous continuum connectivity per- 
colation studies. In what follows we measure the inde- 
pendent connectivity percolation exponents, test the hy- 
per scaling hypothesis, and compare our results to lattice 
connectivity percolation. 

In Fig. [2] we show the data used to determine the 
connectivity percolation exponents in d = 2 for both 
monodisperse and bidisperse [14] systems. Every data 
point is averaged over 5000 trials, each consisting of an 
initially random, dilute state that is subsequently com- 
pressed to the percolation onset at 0p. The correlation 
length exponent v is measured using the standard finite 
size scaling analysis [11], where the root mean square 
width A of the </>p distribution is related to the system 
size L by A ex L x l v . In Fig. [21(a) we plot A 2 versus 
TV and observe a power law behavior with an exponent 
of —2/dv. This gives v = 1.92 ± 0.03, which is larger 
than the value of 4/3 for lattice connectivity percolation. 
Next, we determine the exponent r by directly measur- 
ing the cluster-size distribution. In Fig. EJb) we plot 
n s ((/)p) for a system with 1024 particles and find that 
n s ((/)p) oc s~ T with r = 2.04 ± 0.04. This is consistent 
with the lattice result of 187/91. Finally, in Fig. ^c) 
we plot the average number of particles in the percolat- 
ing cluster iVp er as a function of TV at 0p. We find that 
7V per oc N D / d , where D = 1.88 ± 0.04 is the fractal di- 



present work d = 2 1.92 ±0.03 2.04 ± 0.04 1.88 ± 0.04 

lattice d = 2 [11] 4/3 (1.33) 187/91 (2.05) 91/48 (1.90) 

continuum d = 2 [12] 1.34 ± 0.02 2.0 ± 0.1 1.91 ± 0.04 

present work d = 3 1.10 ±0.03 2.20 ± 0.03 2.56 ± 0.06 

lattice d = 3 [11] 0.88 2.18 2.53 

continuum d = 3 [13] 0.94 ± 0.20 2.22 ± 0.08 2.21 ± 0.56 

TABLE I: Critical exponents for connectivity percolation. 

mension at percolation. This is also consistent with the 
lattice result of 91/48. 

A similar analysis has been applied to d — 3; the mea- 
sured exponents in d = 2 and 3 are collected in Tabled In 
both cases D and r are consistent with lattice and contin- 
uum connectivity percolation, while v is larger than the 
classical value. The hyperscaling relation D(r — 1) = d 
is consistent with our results. 

We have tested the robustness of the measured expo- 
nents by altering key parameters in the force law and by 
varying the preparation algorithm. First, we varied the 
minimum force C over several orders of magnitude from 
10 -4 to 10 _1 and the range of the attraction £ from 
to 4. Second, instead of compressing, we tested for per- 
colation by quenching to zero temperature at constant 
volume. Third, we redefined the connectivity percola- 
tion transition to require percolation in the x—, y— or 
z— directions or all directions at the same time. In all 
of these cases the connectivity percolation exponents re- 
main the same as those listed in Table HI We also expect 
that the measured exponents do not change if the system 
is allowed to coarsen from successive decompressions and 
compressions or other methods. If the system coarsens 
before it percolates, the length scale over which struc- 
ture is affected has a finite extent (equal to the size of 
the largest cluster after coarsening). Near percolation £ 
is much larger than this length scale and coarsening is 
irrelevant. However, we were not able to directly verify 
this conjecture because of insufficient system size. 

Based on these, considerations we conclude that the 
measured exponents (shown in Table [Tj) are universal for 
attractive, frictionless particulate systems. Thus we have 
identified a new universality class — connectivity percola- 
tion with force balance — that is different from connectiv- 
ity percolation in the absence of forces. The universality 
class changes abruptly when Y > (emergence of ex- 
cluded volume) and C > (emergence of attraction). 

Rigidity Percolation The rigidity percolation transi- 
tion at (/>r, which separates the percolated-unjammed 
and jammed states, occurs when a set of particles with no 
floppy modes has percolated. We refer to this set as the 
backbone. In simulations we identify backbone particles 
as those with at least one non-zero contact force (stressed 
particles) . Below (/>r we always observe that no particles 
are stressed and above </>r a percolating cluster of stressed 
particles with no floppy modes exists. Figs. W[b) and (c) 
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TABLE II: Critical exponents for rigidity percolation. 

are screenshots of two systems at rigidity percolation. 
We have also used P > to define the onset of rigidity 
percolation, and obtain identical results [l5j |. 

Similar to connectivity percolation, there are three in- 
dependent critical exponents in rigidity percolation, and 
we choose £, r and (backbone fractal dimension) to 
characterize our system. Numerical values of these expo- 
nents are listed in Table [III along with the corresponding 
lattice results. Rigidity percolation on a lattice has been 
studied in d = 2 [16] and d = 3 0- In the d = 2 lat- 
tice studies the transition is second-order but only the 
exponents and v have been measured. Our measure- 
ment of v differs from the lattice value, while remains 
the same. In the d = 3 lattice studies, the order of the 
transition depends on lattice type, which makes compar- 
ison with our off-lattice results difficult. Hyperscaling is 
consistent with our results both in d = 2 and 3. 

We have also measured the fractal dimension D of 
the entire percolating cluster, and we find that D = d. 
Therefore, while the rigidity percolation transition is sec- 
ond order in the backbone particles, it might appear to be 
a discontinuous first-order transition if the entire perco- 
lating cluster is mistaken for the rigid percolating cluster. 

Extension to finite temperature In this letter, we have 
considered attractive particulate systems at zero temper- 
ature. We are currently investigating the effects of ther- 
mal fluctuations. When T > 0, contact formation due to 
particle or cluster diffusion and contact breaking due to 
thermal activation become relevant. If T is much smaller 
than the energy needed to break a contact and the diffu- 
sion length scale is small, we expect the T = results to 
hold. This has been verified for connectivity percolation 
in our preliminary studies of these systems undergoing 
Brownian dynamics. 

Experimental verification of our predictions can also be 
found in the literature. First, in Ref. [l8j |. experiments 
on colloidal gels demonstrate that the connectivity per- 
colation transition (associated with gelation) occurs at 
a lower density than the rigidity percolation transition 
(associated with arrested dynamics at all length scales). 
This is consistent with our observation that connectivity 
and rigidity percolation occur at different densities. Sec- 
ond, the critical exponents for gelation have been mea- 
sured as v = 1.17 ±0.10 0, r = 2.15 ±0.05 in col- 
loidal suspensions, and v = 1.22 ± 0.16, r = 2.26 ± 0.08 
in silica gels 
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1221 ]. These values are consistent with 
the d = 3 connectivity percolation exponents measured 
here for attractive particulate systems, but are inconsis- 
tent with the connectivity percolation exponents that do 
not include interparticle forces. 
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